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METHOD FOR ESTIMATING CHANGES IN PRODUCT LIFE 
FOR A REDESIGNED PRODUCT 

CROSS REFERENCE TO RELATED APPLICATION 

The present application claims the benefit of U.S. provisional patent application 
Serial No. 60/153,007 entitled "HALTPLUS", filed on September 9, 1999. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention generally relates to reliability engineering. It particularly 
relates to estimating product failure rates using accelerated stress testing data. 

Background Art 

An important design characteristic of any commercial product (particularly 
electronic products) is the reliability of the product (system). The military broadly 
defines reliability as "the probability that an item will perform a required function 
without failure under stated conditions for a stated period of time." (see MIL-HDBK- 
72 1C, "Definition of Terms for Reliability and Maintainability"). Other related 
references on this subject matter include MIL-HDBK-217, "Reliability Prediction of 
Electronic Equipment", MIL-HDBK-883, 'Test Methods and Procedures for 
Microelectronics", and Nelson, "Accelerated Testing: Statistical Models, Test Plans, and 
Data Analyses". 
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Indeed, the commercial viability of the product may be strongly determined by 
the product's reliability as potential users will not spend a significant sum of money on 
an unreliable product. Also, the vital importance of product reliability spurred the 
emergence of Reliability Engineering as an engineering discipline that devotes itself to 
5 establishing, maintaining, and evolving the reliability of a product. 

An exemplary useful reliability parameter is the product's failure rate (Failure In 
Time - FITS, failures per billion device hours) or the inverse measurement of Mean- 
Time-Between-Failures (MTBF). A failure is commonly defined as any event that 
prevents the product from performing its specified operations. A product's failure rate is 

10 the average rate at which the product will fail over its lifetime (e.g., 2 failures over 50 
years). Conversely, the MTBF is total operating time divided by the number of failures 
over that period (e.g., 25 years/failure or 219,000 hours) which is generally regarded as 
the length of time a user may reasonably expect the product to work before a failure 
occurs. Therefore, a primary design objective is to have a very low product failure rate or 

15 inversely a very high MTBF. Indeed, most commercial telephone switching equipment is 
designed for less than two hours of downtime in 40 years, undersea telephone systems are 
designed for less than three failures in 25 years, and personal computers are designed for 
a MTBF of at least 1,000,000 hours (assuming end-of-life replacement of component 
parts). 

20 Therefore, an important reliability engineering procedure is to determine the 

MTBF for a product. Many reliability engineering methods determine MTBF by 
observing the product in the field during normal usage and tracking the frequency of 
failures. However, this method is very costly and time-consuming and does not meet the 
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practical commercial needs of a product manufacturer who wants to accurately predict 
the MTBF early in the development process which is not possible with currently 
available methods. Early prediction of MTBF allows the product manufacturer to make 
any necessary design changes at a low cost to ensure a product deployment with a very 
5 low failure rate. Without early prediction of MTBF, a product manufacturer must rely on 
very time-consuming field use data (requiring months and years of field testing) which 
delays product deployment and adversely affects future sales. Alternatively, product 
quality may be reduced if deployment cannot wait for the field use data or the field use 
data is no longer useful once it is finally obtained thereby increasing the risk of 
10 distributing defective products to the market. 

To provide early prediction of MTBF, product manufacturers commonly use a 
procedure standardized by the military described as the "Bill of Materials" (BOM) 
approach. (Standard MIL 217 and MIL-HNBK-2 1 7F Reliability Prediction of Electronic 
Equipment). The BOM approach theoretically determines the product's MTBF by using 
1 5 the MTBF of each product component. For example, each individual resistor, integrated 
circuit and other component in a system may have a known MTBF, and the MTBF of the 
complete product can be calculated using these values. However, the BOM approach is 
notoriously unreliable because of interactions between product components and other 
factors which may include, but are not limited to manufacturing process control, end-use 
20 environment, and design validation and ruggedness. 

Thus, there exists a need to provide an accurate and early estimation of a 
product's MTBF enabling product manufacturers to deploy relatively defect-free 
products. 
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SUMMARY OF THE INVENTION 

The present invention overcomes the previously mentioned disadvantages by 
providing an efficient method for early estimation of MTBF for a product from 

5 "accelerated stress testing" data. The accelerated stress testing measures the response of 
the product operating in a first environment which is used to estimate the MTBF for the 
product operating in a second environment. In accordance with embodiments of the 
present invention, the first environment may be one or more test environments more 
likely to cause the product to fail than the second environment which is the normal 

10 operating environment for the product. The method disclosed herein advantageously may 
provide an efficient and early estimation of MTBF enabling necessary design changes 
and reduction of development costs. An additional feature of the present invention 
provides an efficient technique for determining the life impact of a product redesign 
wherein the accelerated stress testing data may include data from a plurality of test 

15 environments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram showing an illustrative product development cycle. 

Fig. 2 is a diagram showing linear acceleration factor in accordance with an 
embodiment of the present invention. 
20 Fig. 3 is a diagram showing relationship between accelerated and field use times 

in accordance with an embodiment of the present invention. 

Fig. 4 is a table showing MTBF estimation in accordance with an embodiment of 
the present invention. 

4 
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Fig. 5 is a table showing MTBF estimation in accordance with an embodiment of 
the present invention. 

Fig. 6 is a table showing MTBF estimation in accordance with an embodiment of 
the present invention. 

5 

DETAILED DESCRIPTION 

OVERVIEW 

FIG. 1 provides an illustrative example of a typical product cycle to which the 
present invention may be advantageously applied. It is noted that the product cycle 

10 shown is exemplary and numerous variations thereof may exist dependent on several 
factors, an example being product management procedures. The present invention is 
applicable such variations as well. Thus, it is to be understood that the described product 
cycle is provided for explanatory purposes only and does not serve as a limitation on the 
scope of the present invention. 

15 Typically, a product may be initially conceived and developed during the research 

and development stage. Following a successful proof of concept or feasibility study, the 
product enters the marketing stage where it is commonly marketed to potential customers 
and/or subject to detailed market research assessments. After favorable marketability 
studies are performed, the product may be manufactured and then sold and distributed to 

20 customers in the manufacturing and sales phases. Commonly, reliability testing 

performed during the manufacturing phase uncovers product defects which must be 
corrected. This can result in costly delays to deployment. In embodiments described 
herein, the invention may be advantageously applied prior to the manufacturing stage 
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(e.g., the research and development stage) to discover latent product defects thereby 
enabling efficient and low-cost modification of the product if necessary. 

In accordance with embodiments of the present invention, the method described 
herein uses accelerated stress testing data to predict a Mean-Time-Between-Failures 

5 (MTBF) for a product. An exemplary type of accelerated stress testing that may be used 
is Highly Accelerated Life Testing (HALT). HALT is commonly used in conjunction 
with Highly Accelerated Stress Screen (HASS), a program used during manufacturing 
testing to uncover latent product defects and establish product quality control. HALT 
typically comprises a number of different stress regimens (conditions) designed to 

10 produce failures resulting from product usage in extreme operating environments. HALT 
advantageously provides a measure of the response of the product when subjected to the 
extreme operating conditions. The measured response may advantageously include the 
measured characteristics of the product wherein these characteristics may include, but are 
not limited to gain, operating time, bandwidth, frequency, or other product 

15 characteristics. 

Generally stated, the purpose of accelerated stress testing is to break the product 
so it can then be redesigned and improved. These different stress regiments (tests) can 
include, but are not limited to multiple axis vibrational tests, rapid temperature 
transitions, high/low temperature limits, voltage margining, frequency margining, 

20 electrostatic discharge test (ESD), humidity tests, and any other test designed to find 
weak links in the product design. Commonly, HALT testing is carried out in 
environmental test chambers as part of the effort to uncover product weaknesses and then 
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improve the product. These and other features of HALT and HASS are known to those 
skilled in the art and therefore not described in great detail. 

In an embodiment, an additional feature of the invention provides a method for 
estimating the amount of change in product life resulting from the redesign of a product. 

5 The redesign, for example, may have been spurred by accelerated stress testing data (e.g., 
HALT). This change in product life is referred to as Relative Life ( R). It is to be 
understood that these and other terms used herein are solely used for purposes of clarity 
and should not be construed as any limitation upon the present invention. 

The amount of change in product life can be described either in relative terms or, 

10 if an initial estimate of product MTBF exists, in terms of MTBF improvement. The 
change in product life, R, is relative to some unknown initial product life. This relative 
change represents a proportion (factor) increase or decrease in product life due to 
implemented design changes. The method described herein can also be used to estimate 
the absolute change in MTBF if an initial MTBF estimate is available. An initial MTBF 

15 is that which is characteristic of the product prior to implementation of any design 

changes. Relative Life multiplied by an initial MTBF estimate will yield a new MTBF 
estimate reflecting the results of design changes implemented during HALT. 

Advantageously, the initial MTBF may be derived using BOM X values, wherein 
X represents the failure rate for different product components, in conjunction with 

20 compatible reliability software wherein the X factors in the reliability software database 
are adjusted from their default values so that they reflect realistic levels of actual field 
failure rate for the product components. For example, if the database shows X=70 FITS 
for a crystal oscillator (an exemplary product component), but recent field results indicate 
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that the actual FIT rate has improved by 50%, then the X for the oscillator should be 
adjusted accordingly (e.g., lowered to A,=35 FITS). 

The method described herein is predicated on the idea that accelerated stress 
testing relates to field use conditions and that there is a mathematical relationship 
5 between the accelerated stress test and field use. Specifically, it is assumed that this 
relationship adequately maps from the accelerated stress test time line to the field use 
time line. In an embodiment, the accelerated stress testing data results from linear 
acceleration. 

In accordance with an embodiment, the type of HALT process to which this 
10 method can be advantageously applied involves two separate HALTs (HALT 1, HALT 
2). The first HALT is performed on the original unit (or a previous design of the current 
product) and the second HALT is performed on the redesigned unit that incorporates all 
necessary design changes indicated by the first HALT. Advantageously, HALT 1 and 
HALT 2 provide the same stress tests at the same stress levels, as applied to the original 
15 unit and the redesigned unit, but variations in the stress testing and stress levels can be 
factored into the MTBF estimation as described herein. Using this method of HALT, 
each stress regimen (test) is terminated once the operational limit is found. Furthermore, 
the derivation of this method assumes that the conditions (i.e., stress limits, chamber type, 
test suites, equipment, etc.) under which the second HALT is conducted are the same as 
20 the conditions of the first HALT. This assumption of similarity between HALT 1 and 
HALT 2 conditions is necessary to provide a valid comparison of the results for HALT 1 
and HALT 2 that will yield an unbiased statistical estimate for relative life. 
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However, in practical applications of HALT, it often happens that some stress 
limits in HALT 2 are more extreme than the stress limits in HALT 1 . For example, it 
may happen that the stress limits in Rapid Thermal Cycling (performed as part of HALT) 
are more extreme in HALT 2 due to the improved design. If this situation occurs, the 

5 application of this method may yield an underestimate of the true relative life. In other 
words, this method may yield a biased statistical estimate that is too small. Thus, if 
HALT 2 is run with limits more extreme than HALT 1 then the estimate obtained from 
the invention method should be interpreted as something less than the true amount of 
improvement in relative life. Lacking further information, the magnitude of the bias 

10 (i.e., the difference between the expected value of the estimate and the true value) is not 
known. Further details regarding this issue are described herein. 

DESCRIPTION OF TERMS 

15 Following is a description of the terms that will be used throughout the specification: 

R Relative Life representing the proportion change in product life resulting 

from design changes. R is the parameter of interest (i.e., the true but 
unknown value). 

20 

R The statistic used to estimate R. 

R* The natural log of R. 

25 ~R* The statistic used to estimate R*. 

Ri Estimated change in product life for the i* stress regimen. 



30 



The natural logarithm of Ri . 
MTBF F 2 Mean Time Between Failure on field use time scale for redesigned unit. 
MTBFfi Mean Time Between Failure on field use time scale for original unit. 
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MTBF A 2 Mean Time Between Failure on accelerated time scale for redesigned unit. 

MTBFai Mean Time Between Failure on accelerated time scale for original unit. 

5 tp Failure time on field use time scale. 

t A Failure time on accelerated time scale. 

t A 2i Failure time on the accelerated time scale for the i th stress regimen for the 

1 o redesigned unit. 

t A n Failure time on the accelerated time scale for the i th stress regimen for the 

original unit. 

15 AF Acceleration Factor. 

The standard deviation of j?\ 
at The standard error of R* . 

20 



LINEAR ACCELERATION 

Linear Acceleration, commonly used in single test condition cases, is described 
by the equation tp = AF x t A . In this equation, tp is time on the field use time scale and 
t A is time on the accelerated time scale (see FIGs. 2, 3). AF is the Acceleration Factor 

25 appropriate for the particular stress regimen being used. AF is the time-lengthening 

factor used to determine field use time-to-failure from a time-to-failure measured under 
accelerated stress testing conditions. For example, if the product being tested (e.g., an 
engine) ran 600 hours to failure under accelerated stress testing and has an AF of 3, then 
it is assumed that the engine will operate 600 x3 = 1800 hours to failure under field 

30 operating conditions. It is noted that these acceleration factors need not be known in 
order to apply the method described herein. As will be shown in the following 
derivation, these acceleration factors cancel out of all equations when calculating relative 
life. 
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In applying the method described herein, it is assumed that an unmodified unit 
design (or previous design) can be represented by a population distribution for time to 
failure on the accelerated time scale. Furthermore, a unit with implemented design 
changes can be represented by a potentially different population distribution on the 
5 accelerated time scale. The difference between these two distributions will represent the 
amount of change resulting from the implemented design changes. More specifically, the 
ratio of the observed failure times from each distribution can be used to estimate relative 
change in product life (i.e., Relative Life). 

10 CALCULATING RELATIVE LIFE FOR A SINGLE STRESS REGIMEN 

Relative Life (denoted by R) represents the factor (proportion) change in product 
life resulting from a design change or a group of design changes. As described herein 
for an embodiment, the parameter to be estimated is R. R is defined as the ratio of the 

15 field use MTBF of the redesigned unit to the field use MTBF of the original unit (or a 
previous design). The parameter R is given by equation 1 . 

R = MTBF F2 / MTBF F1 eq. 1 

20 MTBFf2 is the field use mean for the redesigned unit studied in HALT 2 and 

MTBFpi is the field use mean for the previous design studied in HALT 1 . The statistic 
that will be used to estimate this parameter is defined as the ratio of time to first failure of 
the modified design to time to first failure of the previous design. By convention, the 
time to first failure for HALT 2 will always be placed in the numerator and the time to 
25 first failure for HALT 1 will be placed in the denominator. Thus, a relative life value 
greater than 1 will indicate an improvement in product life. A relative life value of less 
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than 1 will indicate a degradation in product life and a relative life equal to 1 will indicate 
no change in product life. It is noted that the calculated statistics for relative life are only 
estimates of the true change in product life and are subject to variability. The method 
described herein also includes a determination of confidence limits for R. 

5 For a specific stress regimen, let tM represent the time to first failure during the 

stress test for the original (unmodified) unit. This value can be advantageously 
determined during the HALT 1 testing procedures. In this notation, the subscript A 
stands for Accelerated and the number 1 indicates the original (previous design) unit (i.e., 
HALT 1). HALT 1 advantageously represents the accelerated stress testing performed on 

10 the original unit (previous design). Let tA2 represent the time to first failure during the 
stress test (i.e., HALT 2) for the modified (redesigned) unit. This value may be 
determined during HALT 2 as HALT 2 advantageously represents the accelerated stress 
testing performed on the modified (redesigned) unit. Similarly, t F i and t F 2 represent the 
times to first failures (hypothetical) under field use conditions for the unmodified and 

1 5 modified units, respectively. 

Estimated Relative Life (denoted by R ) for a single stress regimen can be 
calculated from HALT data using the following equation: 

tA2 

20 R = eq. 2 

tAl 

Using t F 2/tFi as an estimate for R as given in equation 1, equation 2 is derived 
25 from the simple relationship shown in equation 3. 

t F2 AF x t A2 t A2 
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r = = = eq. 3 

t F i AF x t A i t A i 

Relative Life for a given stress regimen as calculated from the above equation 
5 represents an estimate for the relative amount of change in product life resulting from the 
group of design changes implemented as a result of the first HALT. A separate Relative 
Life value can be calculated for each stress regimen (test) in HALT wherein each 
regimen advantageously includes a different type of stress or stimuli being applied to the 
unit under test. An exemplary application of HALT includes five stress regimens (tests), 
10 examples being cold step stress (CSS), hot step stress (HSS), rapid thermal cycling (RT), 
step vibration (VIB), and combined environments (CE). Thus, if there are five stress 
regimens being used in the HALT, there will be five separate values for R . 

Exemplary CSS testing comprises the product beginning operation at a starting 
temperature (e.g., 20 °C ) wherein the temperature is then decreased in successive steps 
15 (e.g., 10 °C steps) until the failure occurs. Similarly, exemplary HSS testing comprises 
the product beginning operation at a starting temperature wherein the temperature is then 
increased in successive steps until the failure occurs. Exemplary RT testing comprises 
applying continuous temperature increases and decreases to the product until failure 
occurs. Exemplary VIB testing comprises attaching the product to a vibration table 
20 wherein the product is subject to vibration in the X, Y, Z and rotational axes. Exemplary 
CE testing applies a pre-determined combination of accelerated stress tests (e.g., RT, 
VIB, etc.) to the product. It is noted there exist several variations of the type of 
accelerated testing used and testing procedures described herein, and so while an 
exhaustive list is not included herein, it is to be understood that any such embodiments 
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using a variation thereof are contemplated here and are within the scope of the present 
invention. 

CALCULATING RELATIVE LIFE ACROSS ALL STRESS REGIMENS 

5 

For calculating relative life, let Ri represent the estimated relative life from the i th 
stress regimen. A composite value for Relative Life (call it R ) can be obtained by 

combining the individual values for Ri from each stress regimen. R will be the estimate 
for R (the true but unknown relative life). Assuming that the time to failure values for 
10 each unit (modified and unmodified units) are exponentially distributed, R can be 

calculated as shown below. In the following notation, k is the number of stress regimens 
in the HALT. 



Ri=^- eq.4 

^Ali 

Let R* =la(R i ) eq. 5 



15 



20 



R =EXP[T]=EXP\-YWKi) eq * 7 
L k m 

As mentioned previously, R is a point estimate for R. 

25 CALCULATING A CONFIDENCE INTERVAL FOR R 

A confidence interval for R can be calculated using the following results, the 

derivation of which is described herein. The standard deviation of R* is given by, 
30 ov = l-81 eq. 8 
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It follows that the standard error of based on k stress regimens is, 




eq. 



It can be shown that the distribution of jf is approximately normal. Let Zo/ 2 
represent the standard normal deviate having a/2 of the area in the tail of the distribution. 
A (l-a)*100% confidence interval for R* is given by, 

zT+~ x^r. eq. 1( 

For example, a 90% confidence interval for R* based on k=5 stress regimens is given by, 



This confidence interval for R* can be converted to a confidence interval for R by 
simply taking the antilog of the limits shown above. Thus, a 90% confidence interval for 
R is given by, 



DERIVATION OF RESULTS 

The following relates the derivation of the results used in the method described 
above. This derivation assumes two independent random variables (X & Y) each of 
which have an Exponential Distribution with parameter X. In order to simplify the 
notation, Y is being used in place of t A 2 and X is being used in place of t A i. 



^ ±(1.645x.81) 
±1.33 




eq. 11 
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The Exponential Density Function is given by, 

f(x)=Ve" X ' x 

The Gamma Density Function defined for x>0 is given by, 

r(x,a,ri- 



»a a - 1 -1 x 
I -x *e 



T(a) 

where, T(a) is the Gamma Function given by, 

r(<x)= 



a - 1 -x, 
x *e dx 



0 

For the special case where cc=l, the Gamma Density Function reduces to the 
10 Exponential Density Function with parameter k. Thus, for a^l , 

T(x, 1 , X)=Exponentiai}^) 

It can be shown through the method of convolution that the density function of the 
ratio of two independent random variables (X & Y) each having a Gamma Density with 
15 parameters oti X and a2,A, respectively, is given by the following density function. Let z 
be a random variable where Z = Y/X. The density function for z is (where 0<z<infinity), 



20 I\a 1 +a 2 -z 

f(3= 



a 2~ 1 



r(a 1 ).r(a 2 V( Z+ i) ai " ha2 



where f(Z) is a Beta Density Function. 

In the special case where Y and X are exponentially distributed both oci=l and 
25 a2=L In this case the density function for the ratio reduces to, 



16 



Docket No. 113337 



f(z)=- 



(z+1) 



From this, it follows that the distribution function for z is, 



F(z)« 



f(z)dz 



(z+1) 



-dz 



1=1 



1 



(1+z) 



10 



15 



L e t w = ln(Z). Using the method of "change of variable", the distribution of W is 
derived as follows: 

G(w) = P(W<w) 
= P(ln(Z)< w) 
= P(Z<e w )-F(e w ) 

Therefore, 

G(w) = F(e w ) = l-(l+e w ) _l 



20 



Given the distribution function G(w), the density function for w can be 
determined as follows: 



g (w)=— G(w) = U+e~; *e = 
dw 



\-2 w 

w\ w e 



l+e v 
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The density function for W can now be used to describe the mean and variance 
for W. Let the mean for W be denoted by E(w) and the variance for W be denoted by 
VAR(w). 



E(w)= 



wg(w) dw 



e exp(w) 
w dw->-ln(l i-exp(w)) + ln(exp(w))'- 



(1 -hexp(w)) 



Evaluating this expression from minus infinity to plus infinity gives the expected 
value for W. 



E(w)= 



wg(w) dw= 



w dw= 0 

d + e w ) 2 



10 Taking the antilog of this result gives an expected value of 1. This is what is 

expected in the special case where both distributions have the same MTBF (i.e., the same 
X). 

15 Similarly, the variance for W is given as follows: 



VAR(w)= 



w 2 - dw -> VAR(w)= " W _ - 2- dilog; 1 + exp( w) ) - 2-li<exp(w))-ln(l + exp(w)) + li<exp(w)) 



l + c ; 



( 1 + exp(w)) 



Where, d ilog(x) = )^dt 
i 1-t 
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Evaluating VAR(w) gives the following, 



VAR(w)= 



(w-E(w)) g(w)dw= 



2 e 



w 



l^e v 



, 1 2 
-dw=-i 

2 3 



Thus, the variance of W is 3.289. Therefore, the standard deviation of W is 1.81 
as stated in equation 8 in the section on calculating a confidence interval for R. This 
result for the variance of W has been shown for the special case where X is the same for 
both Y and X. It can be shown that the same result holds when X is not the same for Y 
10 and X. 

CALCULATING MTBF FOR THE REDESIGNED UNIT 

An MTBF estimate for the redesigned unit can be calculated if an initial MTBF 
estimate for the unmodified unit is available. An initial MTBF for the unmodified unit 
15 can often be obtained via a parts count method (e.g., Bill of Materials - BOM). MTBF 
for the redesigned unit can be calculated as follows: 



MTBF Ri =Rx MTBF 

Redesign initial 



eq. 12 



20 Confidence limits for MTBF can be calculated by substituting the previously 

calculated limits for R into the above equation. 

ESTIMATION OF PRODUCT LIFE 
25 Application of the method described herein to actual HALT results for three 

products (Unit A, Unit B, Unit C) is shown in FIGs. 4-6. As shown in the figures, five 
exemplary HALT stress regimens that were used include CSS, HSS, RT (Rapid 
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Thermal), Vib (Vibration Step Stress), and CE (Rapid Thermal combined with Vibration 
Step Stress). Actual results are shown for three different units and all failure times are 
measured in minutes. The BOM MTBF is used as the MTBF for the original product 
(previous design). 

5 For Unit A (see FIG. 4), the estimate for relative life indicates that the unit life has 

improved by a factor of 1 .57. Generally, it may be estimated that the redesigned unit is 
57% better than the original unit in terms of product life. It is recognized that 1 .57 is 
only an estimate for R and that the true value for the amount of change is most likely 
different from this estimate. The 90% confidence limits give upper and lower bounds on 

10 where the true value for R probably lies. In this example, there is a 90% probability that 
the true value for R is larger than .41 and smaller than 5.93. Since values less than 1 are 
contained inside this interval, it is not possible to conclusively state that the product life 
has improved. 

During the course of accelerated stress testing, it is assumed that only four 
15 physical units are being evaluated during HALT which gives rise to the uncertainty. 
Generally, the failure rate of a product is a function of the product's operating 
environment. As the operating environment becomes harsher, the product's failure rate 
will increase. Therefore, commonly a product's field failure rate is low (<1%) when 
operating in a low-stress environment (e.g., office environment at room temperature) 
20 which consequently makes the probability of finding a failure low as well. HALT testing 
improves the probability of finding a failure without requiring a substantial number of 
test units. Assuming four units being tested, the theoretical probability of finding a 
failure during testing is 25% (one in four units). With experience from previous HALTs 
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and the assumptions described herein, it is determined that this theoretical probability is 
actually higher in practice as the physical sample size of four units actually behaves as 
substantially more than four units when subjected to HALT testing, HALT being an 
exemplary accelerated stress testing method. For example, if a failure mode has a field 
5 failure rate of 1%, then the probability of finding at least one failure in a sample of n=4 
units is 4%. A failure mode refers to the specific type of failure uncovered (e.g., solder 
joint, bad component, etc). If application of HALT testing increases the product failure 
rate (e.g., failures occur more frequently) by a factor of 10, then the HALT failure rate 
2 will be 10% (10 x 1%) and the probability is 34% of finding at least one failure in HALT 
yff 10 using a sample of n=4 units wherein the determination of the increased probability may 
Em be based on probability theory known in the art. At the field failure rate of 1% it would 
ry take a sample size of 42 units to give a probability of 34% of detecting at least 1 failure. 
lS This assumption of the increase in sample size has been generally referred to as the 
!5 "Virtual Sample Size". 

□ 15 Applying the factor 1 .57 to the initial MTBF gives an MTBF for the redesigned 

unit of 468585 hours. Once again, this must be viewed as an estimate. The confidence 
limits for R could be applied to the initial MTBF estimate to give approximate confidence 
limits for the true MTBF for the redesigned unit. 

It is noted that method described herein may be advantageously embodied in a 
20 software computer program compatible with any computer operating system. In an 
exemplary embodiment, the accelerated stress testing data may be used by the program 
(inputted either manually or automatically) to calculate the MTBF (including R) for the 
product. 
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EFFECT OF DISSIMILAR STRESS LIMITS 

It has been asserted that if the stress limits for HALT 2 are more extreme than the 
stress limits for HALT 1 then the estimate for relative life underestimates the true value 
for relative life. The following is a proof of this assertion and is predicated on the 
assumption that more extreme stress limits on a given stress regimen imply a greater 
acceleration factor. Let AF 2 represent the acceleration factor for a given stress regimen 
under HALT 2 and let AFi represent the acceleration factor for the same stress regimen 
under HALT 1. Under the situation where HALT 2 has more extreme stress limits than 
HALT 1, AF 2 is greater than AFi. 

Recall that the ratio t A 2/tAi for a given stress regimen is the statistic used to 
estimate 

MTBF A2 
MTBFm 

It follows, as shown below, that this statistic will underestimate R as given by 
equation 1. 
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Assume AF:> AF 




MTBFfi _ AF2XMTBFA 
MTBFfx ~ ^F: x MTBFa 



MTBFfzxMTBFm AF i 
MTBFfixMTBFa2 AFi 



>1 



MTBFfi ^ 

MTBFfi 



MTBFai 



MTBFa* 



R> 



MTBFaz 

MTBFai 



Since MTBF A 2/MTBF A i is less than the true amount of product life change, the 
statistic used to estimate MTBF A 2/MTBF A i will consistently underestimate the true 
amount of change. 

Although the invention is described herein using primarily HALT stress tests, it 
will be appreciated by those skilled in the art that modifications and changes may be 
made without departing from the spirit and scope of the present invention. As such, the 
method and apparatus described herein may be equally applied to any accelerated stress 
testing designed to expose product defects early in the product development process. 
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What Is Claimed Is: 

1 . A method of estimating a product life, comprising: 

determining accelerated stress testing data for a product, the accelerated stress 
5 testing data representing the response of the product operating in a first environment; 

calculating the mean-time-between-failures for the product operating in a second 
environment based on the accelerated stress testing data. 

2 . The method of claim 1 , wherein: 

said first environment being more likely than the second environment to cause the 
10 product to fail. 

3 . The method of claim 1 , wherein: 

the accelerated stress testing data represents the length of time the product 
operates in the first environment before the product fails. 

4. The method of claim 1 , wherein: 

15 the accelerated stress testing data is derived from a plurality of different stress 

tests. 

5 . The method of claim 4, wherein: 

the different stress tests include a temperature test and a vibrational test. 

6. The method of claim 1 , wherein: 

20 said accelerated stress testing data is derived from either linear, quadratic, or 

exponential acceleration stress testing. 

7. The method of claim 1, further comprising: 

calculating upper and lower confidence limits for the MTBF calculation. 
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8 . The method of claim 1 , wherein: 

said accelerated stress testing data includes bill of materials (BOM) information 
on the product. 

9. The method of claim 1 , wherein: 

5 said step of calculating is performed during the design of the product. 

1 0. The method of claim 1 , wherein: 

said step of calculating is performed prior to manufacturing the product for 
commercial use. 

1 1 . The method of claim 1 , wherein: 

10 said step of calculating is performed using a computer program. 

12. The method of claim 1 , wherein: 

the accelerated stress testing data includes accelerated stress testing data for a 
previous design of the product. 

1 3 . The method of claim 1 2, wherein: 

15 the accelerated stress testing data for the previous design of product is derived 

from stress testing in an environment less likely to cause failure than said first 
environment. 

14. The method of claim 1 2, further comprising: 

calculating a change in MTBF from the previous design of the product. 
20 15. The method of claim 12, wherein: 

said step of calculating includes using the equation EXP [l/k3 k j=i ln(t 2 /ti)]; and 
wherein ti = time to first failure during accelerated stress testing for previous 
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design of the product, and t 2 = time to first failure during accelerated stress testing 
for the product. 

16. The method of claim 12, further comprising: 

calculating a factor increase or decrease in the life of the product as compared to 
5 the life of the previous design of the product. 

1 7. The method of claim 1 2, wherein: 

the accelerated stress testing data is derived from a plurality of different stress 

tests. 

1 8. The method of claim 17, wherein: 

10 the different stress tests include a temperature test and a vibrational test. 

19. The method of claim 12, wherein: 

said step of calculating is performed during the design of the product. 

20. The method of claim 12, wherein: 

said step of calculating is performed prior to manufacturing the product for 
15 commercial use. 

21 . The method of claim 12, wherein: 

said step of calculating is performed using a computer program. 
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ABSTRACT 

A method provides early estimation of product life using accelerated stress testing 
data. In an embodiment, data measured from a product operating in first, high-stress 
environment is used to predict how long the product will operate in a second, normal 
5 operating environment before failure. An additional feature of the present invention 
provides a quantified indication of how much the product has improved from a redesign. 
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Figure 2 - Linear Acceleration 
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Figure 3 - Correlation between Accelerated and Field Use Time Scales 
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FIG. 4 
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FIG. 5 
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FIG. 6 
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IN THE UNITED STATES 
PATENT AND TRADEMARK OFFICE 

Declaration and Power of Attorney 

As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below next to my 

name. 

I believe I am an original, first and joint inventor of the subject matter which is 
claimed and for which a patent is sought on the invention entitled Method For 
Estimating Changes In Product Life For A Redesigned Product, the specification of 
which is attached hereto. 

I hereby state that I have reviewed and understand the contents of the above 
identified specification, including the claims, as amended by an amendment, if any, 
specifically referred to in this oath or declaration. 

I acknowledge the duty to disclose all information known to me which is material 
to patentability as defined in Title 37, Code of Federal Regulations, 1.56. 

I hereby claim foreign priority benefits under Title 35, United States Code, 1 19(a- 
d) or 365(a-b) of PCT or foreign application(s) for patent or inventors' certificate listed 
below or priority benefits under 119(e) of any United States provisional application(s) 
listed below and have also identified below any foreign application for patent or 
inventors' certificate having a filing date before that of the application on which priority 
is claimed: 

U.S. Provisional Application Serial No. 60/153,007 filed September 9, 1999 

I hereby claim the benefit under Title 35, United States Code, 120 of any United 
States application(s) listed below and, insofar as the subject matter of each of the claims 
of this application is not disclosed in the prior United States application in the manner 
provided by the first paragraph of Title 35, United States Code, 1 12, we acknowledge the 
duty to disclose all information known to us to be material to patentability as defined in 
Title 37, Code of Federal Regulations, 1.56 which became available between the filing 
date of the prior application and the national or PCT international filing date of this 
application: 



None 
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I hereby declare that all statements made herein of my own knowledge are true 
and that all statements made on information and belief are believed to be true; and further 
that these statements were made with the knowledge that willful false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of Title 
18 of the United States Code and that such willful false statements may jeopardize the 
validity of the application or any patent issued thereon. 

I hereby appoint the following attorney(s) with full power of substitution and 
revocation, to prosecute said application, to make alterations and amendments therein, to 
receive the patent, and to transact all business in the Patent and Trademark Office 
connected therewith: 



I also appoint Frank Pietrantonio (Reg. No. 32289), Robert L. Hails, Jr. (Reg. No. 
39702), Stuart D. Smolen (Reg. No. 42161), Gary S. Morris (Reg. No. 40735), and Clint 
Gerdine (Reg. No. 41035) of Kenyon & Kenyon as associate attorneys, with full power to 
prosecute said application, to make alterations and amendments therein, and to transact 
all business in the Patent and Trademark Office connected therewith. 

Please address all correspondence to Mr. S. H. Dworetsky, AT&T Corp., P.O. 
Box 4110, Middletown, New Jersey 07748. Telephone calls should be made to Samuel 
H. Dworetsky at 973-360-8120. 



Robert B. Levy 
Benjamin S. Lee 
Susan E. McHale 
Jeffrey M. Navon 



Samuel H. Dworetsky 
Thomas A. Restaino 
Cedric G. DeLaCruz 
Michele L. Conover 
Rohini K. Garg 



Alfred G. Steinmetz 



(Reg. No. 27873) 
(Reg. No. 33444) 
(Reg. No. 36498) 
(Reg. No. 34962) 
(Reg. No. 45272) 
(Reg. No. 28234) 
(Reg. No. 42787) 
(Reg. No. 35948) 
(Reg. No. 32711) 
(Reg. No. 22971) 



Full name of I s joint inventor: Michael K. Brand 



Inventor's signature 





Residence: Renton, King County, Washington 



Citizenship: United States of America 



Post Office Address: 



4408 N.E. 25 m Court 
Renton, Washington 



98059 
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Full name of 2 nd joint inventtfH Harry W. McLe 




Inventor's signature 

Residence: North Bend, King County, Washington 
Citizenship: United States of America 
Post Office Address: 540 S.E. 9 th Street 



ry w. McLean 
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North Bend, Washington 98045 



